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ABSTRACT 



Aims. We study the collision of magnetized irregularities (shells) in relativistic outflows in order to explain the origin of the generic 
phenomenology observed in the non-thermal emission of both blazars and gamma-ray bursts. We focus on the influence of the 
magnetic field on the collision dynamics, and we investigate how the properties of the observed radiation depend on the strength of 
the initial magnetic field and on the initial internal energy density of the flow . 

Methods. The collisions of magnetized shells and the radiation resulting from these collisions are calculated using the ID relativistic 
magnetohydrodynamics code MRGENESIS. The interaction of the shells with the external medium prior to their collision is also 
determined using an exact solver for the corresponding ID relativistic magnetohydrodynamic Riemann problem. In both cases we 
assume that the magnetic field is oriented perpendicular to the flow direction. 

Results. Our simulations show that two magnetization parameters - the ratio of magnetic energy density and thermal energy density, 
Ob, and the ratio of magnetic energy density and mass-energy density, tr - play an important role in the pre-coUision phase, while the 
dynamics of the collision and the properties of the light curves depend mostly on the magnetization parameter cr. Comparing synthetic 
light curves computed from hydrodynamic and magnetohydrodynamic models we find that the assumption commonly made in the 
former models that the magnetization parameter ag is constant and uniform, holds rather well, if ag < 0.01. The interaction of the 
shells with the external medium changes the flow properties at their edges prior to the collision. For sufficiently dense shells moving 
at large Lorentz factors ( i 25) these properties depend only on the magnetization parameter cr. Internal shocks in GRBs may reach 
maximum efficiencies of conversion of kinetic into thermal energy between 6% and 10%, while in case of blazars, the maximum 
efficiencies are ~ 2%. 



Key words. Magnetohydrodynamics(MHD); 
obj ects : general;X-ray s : general 



Radiation mechanisms:non-thermal; Galaxies:jets; Galaxies:BL Lacertae 



1. Introduction 

Relativistic outflows have been observed extensively in blazars, 
a class of active galactic nuclei (AGN) known to show the 
most rapid variability of all AGNs. Their remarkable char- 
acteristic, flares in the X-ray frequency range, usually have 
a duration of the order of one day (Mar aschi et al. 1 1999. 
ITakahashi et al. 1 2000). With the improved sensitivity of 
XMM - Newton, variability on time scales of kilo-seconds has 
been studied in the source Mrk 421 ( jBrinkmann et al. ^2001. 
2003; iRavasio et ain 2004). and recently the spectral evolution 
of this object down to time scales of ^ 100 s could be followed 
dBrinkmann et al. 1 2005). 

Often, the internal shock scenario jRees & Meszaros] l994) 
is invoked to explain the variability of blazars (see e.g., 
[Spadaet al. | 2001; Mimic a 611171 2005, he reafter MAM B05) 
and the early light curves of GRBs jSari & Piranl l995, 
1997; Daigne & Mochkovitch| 1998). One dimensional 
dKino et al. 1200 4; MAMB05) and two dimensional 
dMimica et al. 1 2004, hereafter MAMB04) simulations of 
internal shocks in relativistic outflows performed recently show 
that their evolution is considerably more complex than what can 
be inferred from approximate analytic models: the non-linear 
interaction of two shells leads to a merged shell which is very 



inhomogeneous. Hence, simple one-zone models are of limited 
validity when trying to infer physical properties of the emitting 
region from the flares (resulting from the collisions of shells). 
In MAMB05 we introduced a procedure to analyze spacetime 
properties of the emitting regions in relation to the shape of a 
flare. We showed that under certain conditions one can extract 
flow parameters not directly accessible by cuiTent observations, 
like the ratio of the Lorentz factors of the forward and reverse 
shocks (resulting from the collision of the shells), and the shell 
crossing times of these shocks. 

In the present work we extend our simulations to relativistic 
magnetohydrodynamic (RMHD) flows. To this end we have de- 
veloped a RMHD version of the code presented in MAMB04 
and MAMB05 which we call MRGENESIS. Using this new 
code we performed a systematic study of the influence of the 
initial shell magnetization on the internal shock dynamics and 
on the emitted radiation. Since the new code enables us to treat 
dynamic magnetic fields, we can use it to quantify the accuracy 
of the assumptions made in previous works about the constancy 
and homogeneity of the ratio of the magnetic field energy density 
and the internal energy density in relativistic outflows with inter- 
nal shocks (see e.g.. 



Daigne & Mochkovitch 1998; MAMB04) 
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In Sect.|2]we discuss in detail the initial properties of the col- 
liding shells and describe the numerical method used to simulate 
their evolution. The interaction of two colliding shells and of the 
shells with the external medium are discussed in Sect.[3] The re- 
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suits of our numerical simulations are described in Sect.|4] We 
discuss and summarize our main findings in Sect.|5j 



2. Initial shell properties 

A probable cause of the observed blazar flares is the interaction 
of blobs (or shells) of matter within a relativistic outflow (jet), 
propagating at slightly different velocities. Such an interaction 
happens every time two shells collide after some time depending 
on their relative velocity. Internal shocks produced by the colli- 
sion cause an enhanced emission of radiation which is thought 
to give rise to the observed flares. 

The interaction of the shells is simulated using a rel- 
ativistic magnetohydrodynamic (RMHD) numerical scheme 
jLeismann et al. I 2005) coupled to the non-thermal radia- 
tion scheme developed in MAMB04. The resulting code 
MRGENESIS allows us to simulate the dynamic evolution of 
the magnetic field in a plasma instead of assuming that the field 
is randomly oriented in space and that its energy density is pro- 
portional to the internal energy density of the plasma, as was 
the case in our previous investigation (MAMB05). Thereby we 
are able to check whether the latter assumption, which is widely 
adopted in the literature, actually holds in the course of proto- 
typical two-shell interactions. 

In order to simplify further discussion we introduce two 
magnetization parameters. The ratio of the magnetic energy den- 
sity and the internal energy density of the plasma [[] 
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where F is the bulk Lorentz factor of the flow, and where yad and 
p are the adiabatic index and the thermal pressure of the plasma, 
respectively. 
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is the magnetic pressure in the comoving frame of a fluid element 
moving with a velocity v (and a corresponding Lorentz factor 
F). The magnetic field B is measured in the laboratory (source) 
frame, and we have chosen (here and in the following sections) 
units such that the speed of light c = 1, and that the strength of 
the magnetic field is measured in Gauss. 

We restrict ourselves to the simulation of one dimensional 
models, because we have shown previously that the lateral ex- 
pansion is negligible during the collision of aligned shells, i.e., 
all essential features can be captured using ID simulations 
(MAMB04). We assume that the magnetic field is oriented per- 
pendicular to the direction of the flow. Hence, /Jmag - B^ l{%nT^) 
and cr - B^ liAnpT^). In this special case the comoving mag- 
netic field strength B' is given by the field strength in the 



' Note that there is a typo in the equation defining q-b in MAMB04. 
The left hand side of this equation should read Er with B denoting the 
magnetic field strength in the comoving frame. 
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Table 1. Properties of the simulated models where two shells of 
uniform initial rest mass density p\ - P2 - lO-^^'gcm"-' prop- 
agate through a homogeneous ambient medium (Lorentz factor 
Fext = 2.9, density Pext = lO""' x psheii - lO-^"* g cm-"*) with 
Lorentz factors Fi = 5 and F2 - 7, respectively, ofg and cr" are 
the initial uniform magnetization parameters (ratio of magnetic 
energy density and internal energy density, and ratio of magnetic 
energy density and mass-energy density, respectively), iplp)^ is 
the initial uniform thermal pressure of the plasma in units of its 
rest mass energy density, and B'^^^y^ and B'^^^ give the initial value 
of the comoving magnetic fields of the shells and of the external 
medium, respectively. Note that the magnetic field is oriented 
perpendicular to the flow direction, and that the initial condi- 
tions of model B5 are identical to those of model SIO-FIO of 
MAMB05. 



source frame through the relation B' - B/T. Initially the co- 
moving magnetic field everywhere has the value B' - BIT - 

V4^/(rad-l)aBP = 0.61 ^(a»/10-2)p/(l.lX 10-4 Pext) mG. 

In the latter expression p is assumed to be given in units of PextC^, 
and Pext in turn in units of lO-^-' glcvp?, respectively. 

In all relativistic hydrodynamic (RHD) models simulated 
in our previous studies (MAMB04, MAMB05) the magnetiza- 
tion parameters were (initially) - 10^^ and cr = 10-'', 
respectively. For the RMHD investigation presented here we 
have evolved a set of models whose (hydrodynamic) param- 
eters are similar to those of model SlO-FlO of MAMB05. In 
this (hydrodynamic) model two shells of uniform rest mass den- 
sity 10-^*^ gem--', initial thickness and separation lO''* cm prop- 
agate through a homogeneous ambient medium (Lorentz factor 
Fext - 2.9, density pext = 10-^^* g cm-"*) with Lorentz factors 
Fi = 5 and F2 = 7, respectively. The emitted radiation is com- 
puted using the type-E shock acceleration model with o-e - 10^^, 
- 30, and 7/ = 7 x 10^ (see MAMB04 for details). We sim- 
ulated two groups of RMHD models (see Table 1): 

1. The parameters of models Bl to B6 are identical to those of 
model SlO-FlO of MAMB05 except for the magnetization 
parameter 0;^ which varies from 10-' to 10-^. This subset of 
models allows us to study the influence of the initial magne- 
tization on the light curve and its impact on the dynamics. 

2. In models Tl to T4 the strength of the magnetic field is kept 
constant and equal to that of model B4, but the value of the 
initial thermal pressure of the plasma is varied within a factor 
of ~ 20. 

The relative Lorentz factor between the shells Frei = 
FiFzfl- J(l-rf-)(l -F-2)) fulfifls the condition F,.ei > 



Vl + cr under which of a pair of internal shocks forms (e.g. 
IFan et al. l 2004; |Zhang & Kobayashi| 2005). 
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Although our models are computed in one spatial dimen- 
sion, our simulations are effectively two dimensional. The sec- 
ond dimension is spanned by the energy of the electrons. It is 
discretized into 64 energy bins. As mentioned above, we have as- 
sumed a ratio of shell to ambient density ;i'p - 10^. This value is 
a compromise between those numerically feasible and those ex- 
pected in blazars. Larger density ratios fit the blazar values bet- 
ter, but cause very small time steps. For our value each model re- 
quires about two weeks of computing time on an IBM Power IV 
processor, while a model withxp = 10^ needs about one month 
of computing time. We point out that the value we have cho- 
sen is sufficiently large to account for all qualitative behaviours 
expected for even larger values ofxp- 

3. Interaction of shells with the external medium 

Since the shells move faster than the ambient medium, their front 
edges interact with the ambient plasma. The initially sharp edges 
of the shells evolve into a set of simple waves whose position 
and physical characteristics are calculated by means of an ex- 
act RMHD Riemann solver ( IRomero et al. 2005). For this cal- 
culation we assume that the evolution of the shell edges is one 
dimensional (i.e., along the direction of motion), and that the 
magnetic field is oriented perpendicular to the velocity of the 
shells. According to this set-up, we model blobs hat experience 
a negligible transverse expansion when propagating away from 
the galactic nucleus, and which carry along magnetic fields pref- 
erentially oriented perpendicular to the blob direction of motion. 
A similar situation is encountered when modeling magnetized 
shells emerging from the progenitor system of a GRB except 
that the Lorentz factor of the shells is ~ 100. 

A prototypical evolution of the front edge of a shell is shown 
in Fig.[T] In the particular case displayed here, the external 
medium is assumed to be at rest, while the shell moves with 
F = 17. Both its density and pressure are 1000 times larger than 
in the ambient medium. Inside the shell as well as in the ambient 
medium the magnetic field strength is such that = lO^^*. The 
thermal pressure is 1.11 x 10 "'p everywhere. Under these con- 
ditions two shocks form, a forward and a reverse one, which are 
separated by a contact discontinuity, and which propagate into 
the ambient medium and into the shell, respectively. The struc- 
ture is qualitatively similar to the one found in the pure RHD 
case (see e.g.. Fig. 2 of MAMB04). However, the shocks are 
not purely hydrodynamic ones, but fast hydromagnetic shocks. Q 
Thus, they amplify the magnetic field due to the compression of 
the plasma: the front shock compresses the magnetic field of the 
ambient medium (if non-zero), while the reverse shock amplifies 
the shell's magnetic field. 

The ratio of the magnetic energy density and the thermal 
energy density (thin dotted line in Fig.[T]) decreases in the 
shocked regions. The magnetization parameter cr increases in 
the shocked shell since BjipY) is uniform across shocks and rar- 
efaction when the magnetic field is perpendicular to the veloc- 
ity ( IRomero et al. 2005), i.e., cr oc jp oc p. The particular 
case shown in Fig.[T]represents the qualitative evolution of a sin- 
gle shell. However, as the quantitative properties of the shocked 
region depend on the Lorentz factor of the shell, we have per- 
formed a parameter study to determine these. The four panels 
of Fig.|2] show the variation of the density, the pressure, and 
the magnetization parameter in the shocked regions (mea- 



^ Note that as the magnetic field is assumed to be oriented perpendic- 
ular to the flow direction, fast and slow magnetosonic waves as well as 
Alfven waves propagate with the same velocity. 
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Fig. 1. Flow structure near the front part of a relativistic shell 
moving through an ambient medium (to the right). The initial 
density and pressure of the shell is 1000 times larger than that 
of the ambient medium. The shell propagates with a Lorentz 
factor F - 17, and the initial ratio of pressure and rest-mass 
energy density is 10""*. Two shocks form, a reverse shock (RS) 
propagating into the shell, and a forward shock (FS) propagat- 
ing into the ambient medium. The shocked shell and the ambient 
medium are separated by a contact discontinuity (CD). The var- 
ious lines show (from top to bottom) the rest mass density (thick 
solid line), the thermal pressure (thick dotted line), the ratio of 
thermal pressure and rest mass energy (thin solid line), the mag- 
netic pressure (thin dashed line), and the ratio of magnetic en- 
ergy density and thermal energy density (thin dotted line) using 
arbitrary units. 



sured relative to the corresponding quantities in the unshocked 
regions) with the shell Lorentz factor Several remarkable points 
can be inferred from Fig.|2] 

1 . The ratio of the thermal pressure in the shocked shell region 
and the initial pressure rises with the Lorentz factor. Even for 
relatively small Lorentz factors ( i 5) the thermal pressure 
is 100 times larger than the initial one in the case of an am- 
bient medium at rest. The pressure of the ambient medium 
rises 100 to 10"* times above its initial value. In the case of a 
moving ambient medium and for shell Lorentz factors i 4 a 
rarefaction instead of the reverse shock forms (pL/Psheii < 1; 
flow structure <— HCS — while for larger Lorentz factors 
the situation is as shown in Fig.[T](<— SCS — >). 

2. A comparison of panels (a) and (c), or (b) and (d) shows that 
the properties of the shocked shell region only weakly de- 
pend on the magnetization of the ambient medium provided 
the latter is initially sufficiently small (ffB.ext ^ lO""*) and the 
shell has Fsheii ~ 5, i.e., an initially weakly magnetized am- 
bient medium has no influence on the evolution of the shell. 

3. The rest mass density in the shocked part of the shell (L- 
region in Fig.lTJ is independent of the Lorentz factor for 
r i 5. The comoving magnetic field B' - B/T is indepen- 
dent of the shell Lorentz factor, too, since B' /p = constant 
across shocks for magnetic fields perpendicular to the ve- 
locity field. Since cb oc p^/p, pL - constant, and /jl is 



^ In this notation 7? represents a rarefaction, C a contact discontinuity 
and <S a shock. The arrows indicate the direction of propagation of the 
rarefaction or of the shock with respect to the contact discontinuity. 
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Fig. 2. Variation of the logarithms of pilpskew, Pu/Pext, pilpsixew, PR/Pexi, cb.l, and ffB.R as a function of the shell Lorentz factor, 
Tsheii- Subscripts L and R denote the shocked regions to the left and right of the contact discontinuity, respectively. The density and 
the pressure of the shell are 1000 times larger than in the ambient medium, and Ofsheii = 10"^ (as in the case shown in Fig. [Til. The 
pai-ameters of the ambient medium are: (a) (Text, ofB.ext) = (1, lO"-'), (b) (Text, ffB.ext) - (3, lO"-'), (c) (Text, aB.ext) = (1,0), and (d) 
(Fext, ffB.ext) = (3, 0), respectively. 



monotonically increasing, qtb.l decreases monotonically in 
the L-region. 

The latter fact must be taken into account in the internal 
shock scenario to properly estimate the strength of the magnetic 
field of the faster shell at the moment when the interaction starts 
(i.e., when its front edge starts colliding with the back edge of 
the slower shell), since the magnetic field will be larger, and qtb 
will be smaller than their corresponding initial values. 

If a value Xp > 10^ were used, the leading edges of both 
shells would develop a Riemann fan with the same qualitative 
structure. However, relative to the shell flow the two shocks 
emerging from the leading edges will propagate much more 
slowly, and accordingly the width of the Riemann fan will be 
smaller, too. In the limit Xp ^ °° the width of the fan tends to 
zero, and the shells behave as rigid bodies propagating through 
the vacuum (which is the typical assumption made in analytic 
models of colliding shells). 

Zhang & Kobayashil(2005) have considered in detail the de- 
pendence of various flow parameters on the relative Lorentz fac- 
tor between the shell and the ambient medium, and on the mag- 
netization parameter cr. They also find that for a sufficiently large 
Lorentz factor, appropriately normalized flow parameters are in- 
sensitive to the Lorentz factor, depending only on the magne- 



tization. For example, the dependence of the ratio of magnetic 
pressure and thermal pressure in the shocked region of the shell 
(see panel (e) in Fig. 1 of their paper) is consistent with the trend 
exhibited in Fig.|2] 



4. Numerical simulations 

As in MAMB05, we simulate shell collisions in one spatial di- 
mension only. This approach is justified as all essential features 
of the collision of aligned shells can be captured in one dimen- 
sion, since the lateral expansion is negligible during the colli- 
sion. 

The computational grid consists of 10"* equidistant zones 
covering a physical domain of length 5 x 10'^ cm. A grid re- 
mapping technique (MAMB04) is applied in order to follow the 
evolution of two shells initially separated by lO'"* cm up to dis- 
tances of 10'^ cm with sufficiently good resolution. 

The energy distribution of the non-thermal electrons is sam- 
pled by 64 logarithmically spaced energy bins, and the syn- 
chrotron radiation is computed at 40 frequency values logarith- 
mically spanning the frequency range from 10'^ Hz to 10^° Hz. 
We use the type-E shock acceleration model of MAMB04, with 
some modifications as explained below. 
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Fig. 3. Same as Fig. [T] but for the flow structure near the inter- 
face of two interacting shells (assuming no evolution prior to in- 
teraction). The initial density, pressure and comoving magnetic 
field are assumed to be the same in both shells. The Lorentz fac- 
tor of the faster (left) shell is 7, and that of the slower (right) 
shell is 5. 

The RMHD conservation laws are integrated and the syn- 
chrotron radiation is computed using MRGENESIS, which is an 
extension of the RMHD version of the RHD code GENESIS 
dLeismann et al. 1 2005). For the adiabatic index of the plasma 
a value yad = 4/3 is assumed. The handling of the non-thermal 
particles and the computation of the non-thermal emission fol- 
lows the procedure described in MAMB04, except for the fol- 
lowing two modifications: 

1. Synchrotron radiation: for each grid zone the dynamically 
evolving magnetic field B (measured in the source frame) 
is transformed into the comoving frame defined by the flow 
velocity v (or Lorentz factor F) of the zone in order to 
compute the synchrotron emissivity according to Eq. 14 of 
MAMB04. However, instead of the magnetic field in the co- 
moving frame B' we insert B' sin 6', where 6' is the angle (in 
the comoving frame) between the line of sight towards the 
observer and the magnetic field vector. In the present work 
the flow is assumed to be propagating towards an observer 
and parallel to the line of sight. Since also the magnetic field 
is assumed to be initially perpendicular to the direction of 
the flow, we get B' = B/T and 6' = n/2. 

2. Acceleration of non-thermal particles: we assume that 
non-thermal particles are accelerated in the shocks that 
are detected in the magnetized flow with the standard 
PPM routines implemented in GENESIS ( |Aloy et al. 1 1999). 
However, we replace the jumps in gas pressure by the jumps 
in total pressure (p* = p + Pmug) in the shock detection algo- 
rithm. 

4. 1 . Magnetohydrodynamic evolution 

FigureU shows an example of the analytic solution of the inter- 
action of two shells, assuming they have undergone no evolution 
prior to their collision. This is an approximation of the more real- 
istic situation where the pre-collision evolution changes the con- 
ditions at the edges of the shells that are interacting, which is the 
case of the detailed numerical simulations presented here and in 
MAMB05. Nonetheless, this approximation captures the qual- 
itative structure of the Riemann fan emerging from interacting 



over-dense shells and, particularly the fact that the density jump 
at the contact discontinuity is expected to be small (which can 
also be seen in the results of our numerical simulations; Figs.|4] 
and|5] where the evolution prior to collision was taken into ac- 
count). 

Figures ID and |5] show three snapshots of the magnetohydro- 
dynamic evolution of the 10 models considered in this work. The 
left column displays the initial state, the middle column the prop- 
agation of the internal shocks through the shells after 500 ks, and 
the right column the break out of the shocks from the merged 
shefl after 2000 ks. 

In the middle column the two internal shocks flanking the 
dominant peak inside the shells can be noticed (corresponding 
to the region between RS and FS in Fig.|3ll as well as the front 
structure of the slower (right) shell, which arises from the inter- 
action of the shell with the ambient medium (see also Fig. [Til. 
For models Tl - T4 the thermal pressure within the front struc- 
ture is different, i.e., larger values of Cg (for a fixed magnetic 
field; see Table 1) yield a greater thermal pressure contrast be- 
tween the slow shell and the shocked ambient medium ahead of 
it. Analytically, it is possible to show that the reverse shock prop- 
agates more slowly through the merged shell as ag decreases. 
This trend is exhibited by the models of the B-series (Fig.|4]l, 
where decreases from 10"^ (model Bl) to 10""* (model B6; 
see Table 1). 

For the models of the T-series the differences in the morphol- 
ogy and the dynamics of the merged shell are not as pronounced 
as for the B-models, because the initial magnetization of the for- 
mer models differs only by a fraction a''^(Tl)/a'^(T4-) » 20, 
while a^(Bl)/a''l(B4) k 10^ However, the slowest shell of 
model T4, which has the largest thermal pressure, develops a 
reverse rarefaction instead of a shock as a result of its interac- 
tion with the external medium (see middle panel of Fig. |5]l, al- 
though this difference does not have any noticeable influence on 
the long term evolution of the merged shell. More relevant in the 
T-series is the near independence of the results on the initial ra- 
tio iplpf. Hence, if the initial shells are sufficiently "cold" (i.e., 
{pipf' t 10"^) the exact value of p does not influence the evolu- 
tion. This implies that it may be impossible to infer the (absolute) 
value of the initial shell pressure from observations. 

The final thickness of the merged shells (after 2000 ks) de- 
pends monotonically on ofg, because the total pressure ptot = 
p{l + Q'B/[2(')'ad - 1)]) increases with a^. Model Bl, which has 
the largest initial magnetization of all models (o-g - 0.1), gives 
rise to the thickest merged shell with ~ 2 x 10'"* cm, which is 
of the order of the sum of the initial sizes of the colliding shells. 
Model B6 produces the thinnest merged shell with ~ 7x10'^ cm, 
which is less than the initial size of a single shell. Hence, the final 
size of the merged shells differs little from the sum of the sizes 
of the two initial shells. This can be understood from the fact that 
the final size of the merged shells results from the competition 
of two factors: the expansion triggered when the colliding shells 
are heated up by internal shocks, and the compression produced 
by the collision. 

4.2. Energy conversion efficiency 

In MAMB05 we monitored the efficiency of the conversion of 
bulk kinetic energy into internal energy during the collision of 
the shells. The knowledge of the efficiency is relevant to make 
predictions about the amount of radiated energy in any model in- 
voking internal shocks as the primary source of the non-thermal 
emission of an astrophysical plasma. Here we extend our pre- 
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Fig. 4. Three snapshots from the evolution of models B1-B6. The thick solid line shows the logarithm of the rest-mass density, 
the dashed Une the thermal pressure, the thin sohd hne the comoving magnetic field energy density, and the dash-dotted fine the 
Lorentz factor, respectively. The panels in the first column display the situation 1.5 ks after the start of the simulations, the second 
column after 500 ks, and the third column after 2000 ks (all times are measured in the source frame). The size of the horizontal axis 
corresponds to a length of 4 x lO^'* cm. 




Fig. 5. Same as Fig.|4] but for models T1-T4. The snapshots from models B4 and B5 are shown here again to allow for an easier 
comparison between the models of the T and B series. 
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vious study to magnetized shells, i.e., the bulk kinetic energy 
of the shells can be converted either into thermal energy or into 
magnetic energy. 

We fist define the total instantaneous efficiency of converting 
bulk kinetic energy into magnetic and thermal energy as 



e(T) 1 - 



Jnz,Tmz,T)-l)p(z,T)dz 

j r(z,o)(r(z,o)- i)p(z,o)dz ■ 



(4) 



Next, using the standard relativistic magnetohydrodynamic 
expression for the total energy density measured in the labora- 
tory frame t := p/ir2-;?-;?,™g-pr (where /"i := i+rad7V((rad- 
l)p) + o" is the total enthalpy per unit of mass including the 
magnetic contribution), we can identity the terms correspond- 
ing to the kinetic, the thermal and the magnetic energy densities. 
Taking into account that initially there is not only kinetic energy 
in the flow but also thermal and magnetic energy, we define the 
efficiency of converting kinetic energy into thermal energy as 



[7^d - 1 

Tad 



■P(z,0)- 1 p{z,Q)dz 



I (rad-r 
/ j r(o,r)(r(o,r)- i)p(o,r)dz 



(5) 



and the efficiency of converting kinetic energy into magnetic 
field energy as 



/( 
-/ 
/ 



1 ] 


B\z, T) 


2r2(z, T)j 


An 






2r2(z,o) 


] An 



dz 



dz 



(6) 



/ r(o,r)(r(o,r)- i)p(o,r)dz 



respectively. The quantities T, p, p, and B that appear in the in- 
tegrals are measured either at time f = T or at time t -Q, and at 
position z. The integrals in Eqs.|5] and |6] should extend over the 
whole domain swept up by the shells until the end of the simu- 
lation, but we restrict the domain of evaluation to that of the in- 
stantaneous computational grid. This restriction is justified since 
the dominant contribution to the energy (either kinetic, thermal 
or magnetic) comes from the shells and the shocked ambient 
medium, which are both always covered by the computational 
grid, regardless of the grid-remapping (see above). The contri- 
butions ignored due to grid-remapping are thus negligible re- 
garding the dynamic evolution of the system, but may eff'ect the 
efficiencies as they are of the same order of magnitude ( i 1%). 
Therefore, we re-define the efficiencies given in Eqs|5]and|6]by 
scahng them with e(T), i.e.. 



eth(r) = 



e;(r) 



and 



^mag 



(T) = 



e(T), 



e(T). 



(7) 



(8) 



-mag 



Figs.|6] and |2] show the evolution of eth (upper panel) and 
(lower panel) for models Bl to B6, and models Tl to T4, re- 
spectively. From the upper panel of Fig.|6] one infers that the 
efficiency decreases with at late times for all models. The 
lower panel of Fig.|6]displaying the conversion of kinetic energy 
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Fig. 6. Efficiency of converting bulk kinetic energy into thermal 
(upper panel) and magnetic field energy (lower panel) for models 
of the B-series. Times are measured in the source frame. The 
thick grey line gives an analytic approximation (see AppendixlAl 
for details). 
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Fig. 7. Same as Fig.|6] but for models of the T-series. The graphs 
for models B4 and B5 are added for comparison. 



into magnetic field energy confirms that models with < 0.01 
behave qualitatively differently from models with q-b > 0.01, 
as the global maximum of the efficiency of the latter models is 
shifted to earlier times than in models with ub < 0.01 . The upper 
panel of Fig. |7] shows that the thermal efficiencies eth of models 
Tl to T4 are similar to that of model B4, and that the evolution 
is almost identical for all four models. The efficiency e^ag (lower 
panel) is generally larger for a lower initial thermal pressure (see 
Table, 1). 

In Appendix|A]we give the details of an analytic approxima- 
tion to compute the efficiencies. The approximation takes into 
account that the shells are compressible ^M/c/i and not solids that 
collide inelastically, as has been widely assumed in the literature. 
The thick grey line in Fig. |6] shows our analytic approximations 
for fth (upper panel) and e^ag (lower panel), respectively. The 
parameters for the analytic model are the same as in model Bl. 
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4.3. Light curves 

For each model two light curves were computed, one in the en- 
ergy range 0.1-1 keV (soft band), and second one in the energy 
range 2-10 keV (hard band). The bands are defined to approx- 
imately match those used bv lBrinkmann et ain (2005) when an- 
alyzing their observations. As described in detail in MAMB04, 
light curves are computed by summing the contributions of the 
emissivity from each grid zone in each time step taking into ac- 
count the light travel time to the observer 

Figure|8] shows the normalized soft (upper panels) and hard 
(middle panels) light curves of our 10 models. The difference 
between the light curves of the B-models (left) and T-models 
(right) is obvious, as well as their similarity for all T-models 
and model B4. The instantaneous hardness ratio (HR), defined 
as the ratio of hard and soft counts at a given observational 
time (see bottom panels) rapidly rises during the collision phase, 
and then decreases with the observational time. Hence, at late 
times (t i 15 ks) the detected radiation becomes softer Indeed, 
for models Bl to B3 the light curves and the hardness ratio 
have a multi-peaked structure, a behavior which has been ob- 
served for Mrk421 dBrinkmann etal. . 2005) or PKS2155-304 



(Zhang etal. 2006). In these models having a high magnetiza- 
tion (ttB > 5 X 10"^) the relative amount of radiation emitted in 
both bands changes rapidly with time. Models with > 10"^ 
show a persistent, moderately high hardness ratio (> 0.01) at late 
epochs (f > 15 ks), which correlates with the fact that the RMHD 
evolution of the merged shell produced by initial shells with 
as < 10"^ is characterized by a long -persisting structure radiat- 
ing dominantly in the soft band (see Sect. l5.3b . T-models, all of 
them with low initial magnetic field (equal to that of model B4) 
do not exhibit a variable behavior but the same smooth trends as 
models B5 and B6. 



5. Discussion 

In the following sections we will summarize the main results of 
our numerical investigation. 



5.1. Magnetohydrodynamic evolution 

MAMB05 found that the relativistic hydrodynamic (RHD) evo- 
lution of shell interactions can be divided into three phases: pre- 
collision, collision, and post-collision (rarefaction) phase. The 
evolution of magnetized shells also exhibits these three phases, 
but with some differences. 



5.1.1. Pre-collision phase 

Since the magnetic field inside the shell is perpendicular to the 
direction of motion, the initially sharp forward edge of the shell 
develops into a Riemann fan consisting of three simple waves, 
namely two magnetosonic shocks separated by a contact dis- 
continuity (Fig.[T]i- The reverse shock compresses the shell as it 
propagates towards the shell's rear edge. The compression factor 
(ph/pshew) is almost independent of the Lorentz factor leveling 
off" to a value ~ 10 when Fsheii ~ 5 (Fig.|2|i. As the ram pressure 
exerted by the ambient medium is proportional to pT^, the pres- 
sure increase in the shocked regions rises with the Lorentz fac- 
tor faster than the density yielding a strong heating of the shell's 
edge. Since the ratio of the magnetic field and the rest mass den- 
sity remains constant on both sides of the contact discontinuity, 
the magnetization parameter cb.l is reduced in the shocked am- 
bient medium by up to two orders of magnitude for high Lorentz 



factors (Fsheii i 100; see panels (a) and (c) in Fig.|2]). This de- 
crease of fffi.L is important in those cases where the shells take 
a long time to collide, since the reverse shock (originating from 
the interaction with the ambient medium) can traverse a substan- 
tial part of the faster shell before collision. 

While shell heating also occurs for non-magnetized shells, 
the pre-collision evolution of magnetized shells is further char- 
acterized by an increase of the magnetization parameter ctl in 
the shocked part of the shell (see Sec.O, because ctl is propor- 
tional to the rest-mass density when the field is oriented perpen- 
dicular to the direction of motion. On the other hand, the ratio 
of magnetic energy density and thermal energy density ob.l de- 
creases as the jump in the thermal pressure in the shocked re- 
gion is larger than the jump in density, and thus also larger than 
the jump in the magnetic field. Consequently, at the onset of the 
collision phase the magnetic field of the faster shell is larger 
but dynamically less important than the shell's initial field. The 
properties of the shocked shell (magnetization parameters, den- 
sity contrast with respect to the shell, pressure jump, etc.) de- 
pend only weakly on its Lorentz factor when the latter becomes 
large. This fact is of special relevance in the case of GRB af- 
terglows, which may arise from the evolution of a single shell 
moving initially with a Lorentz factor F ~ 100. The properties 
of the shocked shell are then largely determined by the magne- 
tization parameter cr and its density contrast with respect to the 
ambient medium. 



5.1.2. Collision phase 

Figs.m and |5] (middle column) demonstrate that the ratio of the 
magnetic field energy density (thin solid line) and the pressure 
(thick dashed line) has decreased at the collision phase relative 
to its initial value (first column), consistent with what is seen in 
Fig. [3] (the thin dotted line, denoting a^, decreases in shocked 
regions). The spatial distribution of the thermal pressure in the 
region limited by the internal shocks (created during the colli- 
sion) is very similar for all models. However, the shock structure 
of the slower shell resulting from its interaction with the ambient 
medium is characterized by a larger jump in the thermal pressure 
jump in models with larger initial magnetic fields. 

For the conditions considered in our models, which corre- 
spond to the typical values expected in blazars, maximum effi- 
ciencies of conversion of kinetic into thermal energy of ~ 2% 
are found. Much smaller values are computed for the conversion 
of kinetic into magnetic energy (~ 0.2%). According to our ana- 
lytic model (see Appendix A), the maximum efficiencies of con- 
version of kinetic into thermal energy may rise up to ~ 10% for 
the typical conditions expected for internal collisions in GRBs, 
while the efficiency of conversion of kinetic into magnetic en- 



ergy decreases to fmag i 8 x 10 



5.1.3. Post-collision phase 

Once the internal shocks break out of the merged shell, they 
rapidly stretch it, lowering its mean density and, in general, pro- 
ducing rarefaction waves in the rear part of the merged shell 
(similar to what MAMB05 found in the pure RHD case). At 
source frame times later than x 1500 ks corresponding to a few 
sound crossing times of the merged shells the lateral expansion 
of the shocked regions can be significant, i.e., the synchrotron 
luminosity is overestimated, and the size and thermodynamic 
properties of the merged shells have to be considered with care. 
Indeed, overestimating the pressure and the density in the region 
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Fig. 8. Normalized soft (top) and hard (middle) light curves, and instantaneous hardness ratio (bottom) for models Bl - B6 (left), 
and Tl - T4 and B2 (right), respectively. Note that all hght curves have been normahzed individually. 



bounded by internal shocks could be the reason for some features 
observed at very late times in our synthetic light curves and in 
the hardness ratio (e.g., the secular increase of the hardness ratio 
in the T-models). 

5.2. Evolution of the magnetic field 

The evolution of the magnetization parameter cb at the point 
of maximum magnetic field at any given instant of time in the 
source frame shows two epochs separated by a sharp disconti- 
nuity at about 300 ks (Fig.|9]l. The instantaneous location of the 
maximum magnetic field is indicative of the region which instan- 
taneously (in the source frame) contributes most to the radiated 
energy. Hence, it is a representative point of the binary shell evo- 
lution. Left of the discontinuity the maximum magnetic field is 
produced by the reverse shock (caused by the interaction with 
the ambient medium) of the faster shell prior to the collision as 
it compresses the fluid in the shell. The part right of the discon- 
tinuity corresponds to a point inside the reverse internal shock 
formed during the collision phase (see Sec. 15. II ). The sharp dis- 
continuity in QfB occurs when the global maximum of the mag- 
netic field shifts from the external shock region of the faster shell 
to the internal shock region of the merged shell. 

All curves shown in the upper panel of Fig.|9] are qualita- 
tively very similar, especially regarding the relative magnitude 
of the sharp drop in at 300 ks, which is almost identical for 
all the B-models. This implies that the system does not com- 
pletely "forget" its initial conditions. Although inside the region 
shocked by internal shocks is about one order of magnitude 
smaller than the initial magnetization in the shell, the order- 
ing of the models according to ofg is preserved, at least until 
reaches its maximum value at f ~ 900 ks. 

The second magnetization parameter cr" is non degenerate 
in the B-models (Fig.[TOl upper panel), since during the collision 
kinetic energy is dissipated, increasing the magnetic field in the 




500 1000 

Tliis] 



Fig. 9. Ratio of the comoving magnetic energy density and the 
thermal energy density as a function of time in the source 
frame for the grid point with the maximum magnetic field. The 
upper panel shows the B-models which all have the same initial 
pressure, and the lower panel displays the T-models which all 
have the same initial magnetic field. 



shocked region, while the initially small thermal energy compo- 
nent has very little influence. 

A qualitatively different behavior is observed for the mod- 
els of the T-series (Fig.|9] lower panel). These have a different 
initial magnetization parameter ffg, but the absolute value of the 
magnetic field is initially the same in these models. Thus, qb 
evolves very similarly in these models after the discontinuity in 
ffB is encountered. The evolution of cr is degenerate (see lower 
panel of Fig.fTOli. since all models have the same initial value cr" 
and undergo an almost indistinguishable evolution. 
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Fig. 10. Same as Fig. |9] but showing the time evolution of the 
ratio of the comoving magnetic energy density and the kinetic 
energy density. 



We therefore conclude that, in order to describe the essentials 
of the magnetohydrodynamic properties of binary collisions of 
cold shells, one can safely exclude the initial thermal pressure 
(or, equivalently, the internal energy) as a relevant parameter, 
and may use some canonical value instead, which is sufficiently 
small compared to the kinetic energy. 

5.3. Radiation 

The light curves of all models having the same initial magnetic 
field look similar (Fig.[8j right panels), which is consistent with 
the picture described in Sect. 15. 2[ as during the collision phase 
the initial value of cr is the relevant parameter to account for the 
magnetic field evolution. For large values of cr the hght curves 
are shorter and multi-peaked, while for decreasing cr° they be- 
come single-peaked and have a longer decreasing tail which is 
progressively softer (Fig.[8j left panels). The border separating 
cases with a multi-peaked or mono-peaked light curve is given 
by the combination of the value of cr and a critical value of 
of the order of 10"^ (see Sect.|43Tl. 

Comparing the light curves of the current work with those 
given in Fig. 4 of MAMB05, we find that also in the magnetized 
case hard light curves are shorter in duration than soft ones, and 
that only low magnetization (ai < 10"^) models have a shape 



similar to those of MAMB05, which is not surprising since in 
the limit of Cg — > the pure hydrodynamic limit must be recov- 
ered. As can be seen when comparing the lower panels of Fig. [8] 
with the lower panel of Fig. 4 of MAMB05, the initial rise of 
the hard light curves is not smooth, but exhibits a small 'kink' 
both for non-magnetized and magnetized models. This similar- 
ity indicates that the assumptions made in MAMB05 are valid 
for small cKg, but need not necessarily hold for high initial a^, 
too. 

The maxima of the soft and of the hard light curves, as well 
as the total hardness ratio (ratio of the time-integrated number 
of counts in the hard and in the soft bands) depend on the ini- 
tial magnetic field for the B-models (Fig.fTTI). Our models show 
that the maxima of the light curves have a logarithmic-like de- 
pendence on the initial magnetic field energy density. In order to 
explain the weaker dependence of the light curve maxima and 
hardness ratio on the magnetic field for increasing magnetiza- 



• Bl-B6soft 
□ Bl-B6hard 
O HR total 



log b [arbitrary units] 

Fig. 11. Maxima of soft (filled circles) and hard (squares) light 
curves for models Bl to B6 as a function of the square of the 
initial comoving magnetic field. The total hardness ratio (HR) is 
defined as the ratio of the time-integrated number of counts in 
the hard and in the soft bands, respectively. It is displayed with 
the diamond symbols. 



tion, we approximate the typical high-energy synchrotron spec- 
trum (Fig.[T2]l as a broad peak with a power-law decay towards 
higher photon energies. The position of the peak shifts to higher 
energies for larger values of the magnetic field. Q. For models 
with lower magnetic fields the part of the spectrum in the hard 
observational band is the power-law tail resulting in a strong de- 
pendence of the light curve on the value of the magnetic field 
during the evolution. For the high magnetic field models, the 
maximum of the spectrum partly extends into the hard band, i.e., 
increasing the magnetic field does not change the observed en- 
ergy in the hard band as dramatically as in the previous case. 
This can explain the saturation of the dependence of the hard- 
ness ratio on the magnetic field. 

Here we stress once more that the light curve at any given 
time is the result of the sum of the emission from different points 
radiating at different times in the source frame. This explains 
why, e.g., the spectrum of model Bl (Fig.fTSI) has a more com- 
plicated shape than that of model B6. Nevertheless, our basic 
conclusion about the dependency on the magnetic field should 
hold in both cases. 

A similar study as for the B-models has been done for the 
light curve maxima of models Tl to T4 as a function of the 
initial thermal pressure. A very weak dependence of the radi- 
ation on the initial (small) pressure was found, which is consis- 



For Q-B > 0.01 the spectra of the models become more flattened in 
the soft band, and the global maximum shifts towards lower energies. 
This effect is partly due to the faster cooling of the particles, which ra- 
diate at lower energies as time progresses, and partly due to the fact 
that the spectra emitted by particle populations at later source times 
(having less energetic distributions) are still visible in the observational 
band (due to a sufficiently high magnetic field). This complicates the ob- 
served spectrum in the soft band. Nevertheless, the high energy part of 
the spectrum has basically the same shape as that of the Iow-q-b models. 
Ghirlanda, Ghisellin i & Celotti| (2004) found that a similar mechanism 
(decrease of the peak spectral energy plus an steeper increase of the 
vFy spectrum below the peak) may account for the larger hardness ratio 
of short GRBs. In our case, this behavior is triggered by the magnetic 
field while in lJanka et al. 1 (2006) the effect was attributed to the Lorentz 
factor. 
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Fig. 12. Spectra at the soft light curve maximum for models of 
the B-series. In order to provide an absolute luminosity scale 
we assume a shell radius = 10'^ cm. The edges of the soft 
and hard observation bands are indicated by two pairs of vertical 
lines. 

tent with the discussion of Sec. l5.2l about the importance of the 
initial value of cr and the unimportance of for the collision 
dynamics and the emitted radiation. 

5.4. On the energy subtraction mechanism 

The mechanism used by us to extract energy from the thermal 
plasma and to transfer it to the non-thermal population of radi- 
ating electrons is based on the type-E macroscopic recipe pro- 
posed by MAMB04. In the type-E model one subtracts from the 
thermal plasma a fraction se of the internal energy dissipated 
in shocks. This fraction se is a parameter of our models which 
is poorly constrained by the current knowledge of the physics 
of particle acceleration at shock fronts, a typical value being 
ae ^ 10"^. In this subsection we critically discuss the impact of 
different choices of ae on the result. To this end we computed two 
additional models having the same model parameters as model 
B5, but assuming different values for se: In model B51 we set 
se = 5 X 10"^, and in model B52 we choose se = 10"^ (see 
Fig.inil. 

Regarding the physical parameters of shocked regions (i.e., 
those regions affected by the energy subtraction mechanism) a 
comparison of the evolution of models B5, B51 and B52 pro- 
vides the following insights: 

1 . Increasing the value of se yields smaller values of the thermal 
pressure, i.e., the shocked regions become cooler. 

2. The magnetic energy density grows as se is increased in the 
region swept up by internal shocks. This growth is consistent 
with the fact that a more efficient cooling needs higher mag- 
netic field strengths (because that reduces the cooling time). 

3. The size of the region containing the shocked ambient 
medium ahead of the slower shell (R-region in Fig.[rii in- 
creases as se decreases, because the amount of cooling is 
reduced (point 1 above), i.e., the forward shock is not as 
severely weakened as in models with larger se. Indeed for 
se > 0.01 the forward shock is almost suppressed, and thus 
almost no R-region forms. The increase in size of the R- 
region leads to a decrease of the density, and consequently 
of the magnetic field energy density in that region (Fig. [13] 



panels in middle row). Note that this trend is opposite to the 
increase of the magnetic field energy density observed in the 
internal shock region. 

4. The state reached by models with different se at late epochs is 
quantitatively very different (Fig.[T3]panels in bottom row). 
The energy subtracted from the merged shell (which is swept 
up completely by shocks) is higher for increasing values of 
se. Hence, the merged shells of models B5 and B5 1 are much 
cooler (about one order of magnitude) than that of model 
B52, i.e., the most evident effect of the shell cooling is the 
much smaller expansion of the merged shell at late times. 

5. The larger the value of se, the larger is the magnetization 
parameter a^. In the collision phase, 500 ks after the start 
of simulation (Fig. [13] panels in middle row), the value of 
(Tb at the point having the maximum density (in the internal 
shock region) is more than one order of magnitude higher for 
modelBSl {a^ = 3.4x10"^) than for model B52 (^b = 1.5x 
10""^). The increase of results from the combined effect of 
the decrease of the thermal pressure (point 1 above) and the 
increase of the magnetic field strength (point 2 above). This 
result is consistent with the physical expectation that larger 
magnetic fields (which, in our case, also imply larger values 
of a^) will cool down the shocked regions faster, i.e., for the 
same elapsed time, more energy will be radiated away for 
larger magnetic fields than for smaller ones. 

6. The value of cr at the point of maximum density is rather 
insensitive to the choice of se (for models B5 1 and B52 the 
value of cr at 500 ks is 2.7 x lO"*" and 2.2 x 10"'', respec- 
tively). This strengthens our argument that cr is the appro- 
priate parameter to trace the RMHD evolution of colliding 
shells, as the ordering of the models according to cr remains 
unchanged even if one varies the value of se. 

5.5. Comparison with biazar properties 
5.5.1. Emission properties 

For models Bl - B6 the radiative luminosity of a single flare 
lies in the range lO'*^ - lO'*'* erg s"' assuming a shell radius |3 
of 10'^ cm (Fig.[T2li. This luminosity roughly agrees with the 
lower limit observed for the moderately luminous blazars whose 
synchrotron peak lies in the X-ray band (see e.g.. lFossati et al. I 
1998). As the absolute value of the luminosity depends non- 
linearly on the parameter se (see § I5.4l i. the efficiency of the 
conversion of kinetic energy into radiation is somewhat uncon- 
strained. For a homogeneous cylindrical shell with a rest mass 
density p, Lorentz factor F, and velocity v the kinetic luminosity 
can be defined as 

L^,:^7TRl{hT-\)pTv. 

For our initial models the kinetic luminosity is ^ 2 x 10^^ erg s ' , 
implying a conversion efficiency of the order of 0.01 % -0. 1 % for 
se - 10"^. We emphasize that although both the kinetic and the 
radiative luminosities are unconstrained by our one-dimensional 
simulations up to a factor nR^, their ratio and hence the efficiency 
of conversion of kinetic into radiative luminosities is indepen- 
dent of this factor. 

Our models are set up to have a uniform magnetization cr" 
both in the shells and in the ambient medium. Consequently, 
once the density of the ambient medium is fixed, the magnetic 



' Since our numerical study is one-dimensional, the luminosity of the 
models scales as ttRi, where is the shell radius. 
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Fig. 13. Same as Fig.|4]but for model B5 (middle column) and two additional models: model B51 with ae - 0.05 (left column), 
model B52 with se = 0.001 (right column). 



and 



field strength in the shocked regioij^ roughly corresponds to the 
values inferred from observations of blazars, i.e., ^ 0.01 - 1 G 
(e.g.. lGhiselhni et al. 1 1998). Increasing the ratio of shell to am- 
bient medium density ;i'p, will increase the kinetic luminosity of 
the shells proportionally. Assuming a uniform initial magnetiza- 
tion, an increase ofxp will, however, drive the value of the mag- 
netic field in the shocked regions above the aforementioned val- 
ues. Thus, in order to set up future simulations that mimic more 
powerful blazars with realistic magnetic field strengths, we will 
have to consider denser shells whose magnetization parameter 
cr" is smaller than that of the ambient fluid. 

The number density of radiating non-thermal electrons is 
i 10"* - 10^ cm"^, which is one or two orders of magnitude 
higher than the usually inferred values i 10^* - 10''cm"^ (e.g. 
IGhisellini et al. 1 1998). However, these values are obtained un- 
der the assumption that the emitting region is a sphere, whereas 
in our case the emitting region is a cylinder whose height L ~ 
lO'^^cm is much smaller than its radius /?s ~ 10'^ cm. Hence, 
for the same number of emitting electrons, the number density 
of radiating particles in this smaller volume is expected to be 
larger 




Fig. 14. Normalized cross-correlation function (see Eq.|9l) for 
models Bl to B6. Bins of 16.6 s were used when computing the 
cross-correlation functions. 



5.5.2. Light curves 

Following IBrinkmann et al. 1 (2005). we calculate the cross- 
correlation function (CCF) of the light curves of models Bl to 
B6. If the observed counts are binned into equidistant time 
intervals At, the CCF is defined as follows, 

Ii/(^soft(i AO - Xsoft)(Xhard((' + k)At) 

CCF(kAt) := , (9) 



* The absolute scale of B depends on pex, (Sec.|2ll. Since the value 
of the magnetic field strength in the shocked region (Bshock) may be 
computed either analytically or numerically, we can tune pext such that 
Sshock ^ 0. 1 - 1 G. 



where k - 0, 1, (A^ - 1), and iisoft, and Xhard are the number 
of counts averaged over the N-intervals for photons detected in 
the soft and in the hard band, respectively, o-^of^ and o-jjari are 
the standard deviations of the samples with respect to the corre- 
sponding average number of counts in the soft and hard bands, 
respectively. n{k) is the number of pairs (/, i + k) where Xsoft(iAf) 
and Ji:haid((' + k)At) are both nonzero for a given k. The time in- 
terval kAt is called the time lag. Significant correlation at e.g., 
negative lags implies that soft band variations are observed later 
than hard band variations. 

Contrary to the observations shown by ' Brinkmann et al. I 
(2005; Fig. 3), the maxima of the CCFs are not located at zero 
time lag, but at negative time lags (Fig.O, i.e., our simula- 
tions only predict the existence of soft lags. However, the over- 
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all shapes of the CCF curves are well covered by our models, 
including both symmetric and asymmetric ones. 

For small initial magnetic fields, the CCFs become more 
asvmmetric. lRavasio et al. 1 (2004) pointed out that the asymme- 
try can be reproduced by modeling light curves as initially rising 
linearly, and subsequently decreasing with different time scales 
at different frequencies. This behavior can be understood from 
our models as well. For lower magnetic fields the light curve 
is visible for a shorter time in the hard band than in the soft 
band (Fig.[T2b. thus giving rise to different time scales. We point 
out that these time scales are not only related to the cooling of 
the relativistic electrons, but to the radial expansion of the emis- 
sion region as well, whereby the magnetic field inside this re- 
gion is decreased. The fact that we only see soft lags is probably 
a consequence of the particular acceleration mechanism model 
we have chosen (fixed lower and upper limits of the non-thermal 
electron energy distribution; type-E model of MAMB04). The 
influence of the acceleration mechanism will be the subject of 
future work. 



5.6. Summary 

A detailed study of the interaction of magnetized shells in rel- 
ativistic outflows has been performed. An idealized interaction 
of a shell with an ambient medium was studied using an ex- 
act Riemann solver ( Romero et al. 1 2005). while the full inter- 
action of two colliding shells was simulated in one spatial di- 
mension using a RMHD code (' Leismann et al. ' 2005) coupled 
to the synchrotron emission scheme of MAMB04. The models 
are parametrized to address blazar jets under the hypothesis that 
blazar light curves result from internal shocks within relativistic 
jets. However, some results can be extrapolated to the dynam- 
ics of internal shocks in GRBs, when taking into account that in 
the latter case the Lorentz factors are one order of magnitude or 
more larger than in blazar jets, and that therefore much smaller 
relative velocities are encountered between the colliding shells. 
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Fig. A.l. A schematic representation of the ideaUzed interaction 
of two colliding shells consisting of four distinct regions. The 
shells are assumed to move towards the right. 



1. unshocked faster shell (between xi^(T) and xrs(T)): ph, Pl, 
Bl, Fl 

2. shocked faster shell (between xi^(RS ) and xcoiT))'. Pls, Pls, 
Bls, Fcd 

3. shocked slower shell (between xcd{T) and xps(T)y. pRs, Prs, 
Brs, Fcd 

4. unshocked slower shell (between Xps(T) and xp^(T)): pr, /?r, 
Br, Fr 

The five (time-dependent) boundaries of the four regions are 
from left to right: 



--2 
RS 



-2 
CD 



1. left edge of the faster shell at xi^(T) 

2. reverse shock at XRs(r) = Axl + T - Fj 

3. contact discontinuity at xcoiT) - Axl -i- T ^1 - F, 

4. forward shock at xpsiT) - Axl -i- T ^1 - Fp^ 

5. right edge of the slower shell at xr(T) - Axl + Axr + 

where Axl and Axr are the initial thickness of the faster and 
of the slower shell, respectively. The reverse shock, the contact 
discontinuity and the forward shock originate at the interface be- 
tween the two shells. The origin of the coordinate system coin- 
cides with the left edge of the faster shell at T = 0. 

Knowing the initial values of p, p, B, and F in the shells, we 
use the exact Riemann solver ( IRomero et al. 1 2005) to compute 
values in the intermediate states. We then define the total kinetic, 
thermal and magnetic energy at a time T as follows: 



Appendix A: Analytic modeling of energy 
conversion efficiencies 

In this section we idealize the interaction of two colliding shells 
by neglecting the pre-collision phase, and assuming that the 
dynamics of the collision can be approximated by solving a 
Riemann problem at the interface of the two shells. Fig. lA.ll 
shows an example of the flow structure in this case. From left 
to right we identify the following four distinct regions Q, each 
with its rest-mass density p, thermal pressure p, magnetic field 
B (assumed perpendicular to the direction of propagation), and 
the fluid Lorentz factor F: 



' The fluid in the faster and slower shocked shell regions has the same 
Lorentz factor as the contact discontinuity, i.e., Fcd- 



Ek(T) := FL(FL-l)pL(xRs(r)-XL(r)) 

+Fcd(Fcd - 1)Pls {xcd{T) - XRs(r)) 
+Fcd(Fcd - 1)Prs {xp^ij) - xcd{T)) ' 
+F«(F«-l)pR (xR(r) - XFs(r)) 



£th(r) 



VTad-l / 

+ (^^F^p - l) ;,LS ixcD(T) - XRs(r)) 

-rcD - ij/'RS (■^FS(7') - ■^Cd(7')) 
-F^ -l)/7R (xR(r) - XFs(r)) 



(A.1) 



Tad 

I Tad - l' 
Tad 

I Tad - 1 



, (A.2) 
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and 
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^RS 


CD 




4-K 



(xRs(r) - xUT)) 



(xFs(r) - xcD(r)) 



(A.3) 



As in Eq.m we define the total efficiency for conversion of 
kinetic into thermal and magnetic energy 



e™(r) := 1 



(A.4) 



£k(0) 

the efficiency of conversion of kinetic into thermal energy 

£k(0) ' ^ 
and the efficiency of conversion of kinetic into magnetic energy 



e^"" (T) - 



E mag ( ^) ^ mag (0) 

Em 



(A.6) 



The gray thick line in Fig.|6]shows e™ (upper panel) and e^^g 
(lower panel) for the initial shell set-up of model Bl. We as- 
sume that the shells collide after a time of 335 ks, determined by 
their initial separation (10'"* cm) and Lorentz factors (7 and 5, re- 
spectively). We compute the efficiencies until the forward shock 
reaches the edge of the slower shell, as the analytic approach 
is no longer valid later on. There is a good agreement between 
the final value of e™ (0.0191) and the maximum value of €± for 
model Bl (0.0195). The slope of the initial efficiency increase 
is rather well captured, too. The instantaneous discrepancies be- 
tween the numeric and analytic estimates mostly arise from the 
pre-collision evolution of the shells. 

We have also applied the above method to the conditions ex- 
pected to be found in GRBs. Defining the relative Lorentz factor 
of the shells as Ag = (Fl - rR)/rR, we show e^™ and e,™g in 

Fig. lA. 21 assuming a low magnetization parameter cr = 3 x 10"^ 
and initially cold shells with pip - 10"^. Different curves corre- 
spond to different values of (Fr, Ag). Our results show that, for 
the typical conditions considered for internal shocks in GRBs 
(e.g., |Daigne & Mochk ovitch* 1998), maximum efficiencies be- 
tween 6% and 10% can be obtained. 

Figure lA3] shows the maximum efficiencies (corresponding 
to the final states in Fig. |A.2t for a fixed Lorentz factor Fr as a 
function of the relative Lorentz factor keeping all other pa- 
rameters the same as in Fig. lA.2l The efficiency depends strongly 
on Ag, and is independent of Fr for large values of Ag. 

We have also investigated the dependency of the efficiency 
of conversion of kinetic energy into thermal and into magnetic 
energy on the magnetization parameter cr for two selected sets 
of the parameters (Fr, Ag) which are representative for blazars 
(Fr, Ag) = (5, 2) and GRBs (Fr, Ag) = (100, 1), respectively. As 
we see in Fig. |A.3| increasing the magnetization tends to increase 
at all times independently of the value of the shells' Lorentz 
factor (i.e., both for blazars and for GRBs). Indeed, e™„ reaches 
values of 0.2 - 0.3 for magnetically dominated shells (cr = 10). 
More moderate magnetizations {cr - 1) but still much larger 
than the ones we have considered for the numerical models in 
this work yield e^^g ~ 0.1. The trends for e™ are not monotonic. 



(100, 1) 

-- (100,4) 
(100, 9) 
— (200, 1 
- (200;4 
(200, 9) 




0.4 0.6 
T [normalized units] 

Fig. A.2. Efficiency of converting bulk kinetic energy into ther- 
mal (upper panel) and magnetic field energy (lower panel) as a 
function of time. Shown are six models with different values of 
(Fr,A^) and with cr = 3 x 10"^ where Ag = (Fl - rR)/rR is 
the relative Lorentz factor, and where Fl and Fr are the Lorentz 
factors of the faster and the slower shell, respectively. 




Fig. A.3. Efficiency of converting bulk kinetic energy into ther- 
mal (upper panel) and magnetic field energy (lower panel) as 
a function of the relative Lorentz factor A^, where A^ = (Fl - 
Fr)/Fr, and where Fl and Fr are the Lorentz factors of the faster 
and the slower shell, respectively . Shown are five models with 
diff^erent Fr (5, 10, 50, 100, 200) and with cr = 3 x 10"^ 



and depend on the chosen scenario. For blazar conditions, e^™ 
decreases with increasing values of cr. For the conditions met in 
GRB flows, there is an asymptotic trend to increase e™ for very 
large values of the magnetization parameter (cr > 1). When cr is 
increased, the efficiency of conversion of kinetic into magnetic 
energy becomes larger than the efficiency of converting kinetic 
into thermal energy. We point out that values of cr 1 imply 
that initially the shells may have a magnetic energy larger than 
the thermal energy. When the initial internal energy is larger than 
the magnetic energy, then > e^^g holds , while in the opposite 
case > e^^. 
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